A generic problem with purely metric formulations of MOND 
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We give a simple argument to show that no purely metric-based, relativistic formulation of Mil- 
grom's Modified Newtonian Dynamics (MOND) whose energy functional is stable (in the sense of 
being quadratic in perturbations) can be consistent with the observed amount of gravitational lens- 
ing from galaxies. An important part of the argument is the fact that reproducing the MOND force 
law requires any completely stable, metric-based theory of gravity to become conformally invariant 
in the weak field limit. We discuss the prospects for a formulation with a very weak instability. 
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1. Introduction: A large spiral galaxy might contain 
a mass of M ~ lO^^Af© ~ 10'*^ kg in the form of stars 
and gas. Almost all of this mass lies inside a radius of 
i? ~ 10 kpc ~ 10^''^^^ m. With numbers on these scales 
it is easy to see that the gravitational acceleration ought 
to be Newtonian, 
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Hence neutral Hydrogen in a circular orbit of radius r > 
R should be observed to move more slowly as r increases, 
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As is well known, it does not. In the scores of galaxies for 
which high quality rotation curves can be obtained 0, 
one actually finds that the asymptotic speed approaches 
a constant which is proportional to the fourth root of 
the total luminosity. This is known as the TuUy-Fishcr 
relation jij. 

Dark matter is the conventional explanation for the 
observed failure of the Newtonian prediction jSJl. It was 
originally invoked to provide the mass necessary for New- 
tonian gravitation to bind galaxies to their host clusters 

Q . To explain flat rotation curves it is supposed that 
galaxies are surrounded by halos of weakly-interacting, 
non-relativistic particles whose contribution to the total 
mass density falls off much more slowly than the mass 
density in stars and gas . 

Dark matter has been impressively successful in recon- 
ciling the observed universe with general relativity and 
its Newtonian limit 6]. However, it is not without prob- 
lems. First, there has been no direct detection of the par- 
ticles which comprise it, except for neutrinos which can 
provide only a small fraction of the necessary mass and 
which are too light to cluster on galaxy scales. Second, 
what accounts for galactic rotation curves is an isother- 
mal halo piso ~ l/r^ of dark matter. But numerical 
simulations of dark matter consistently produce the den- 
sity profile of Navarro, Frenk and White 0, which falls 
like pnfw ~ ^/r^ at large distances. Dark matter does 
not explain the TuUy-Fisher relation. Nor can it eluci- 
date the curious fact that the breakdown of Newtonian 



gravity, sourced by only the mass in stars and gas, seems 
to occur at the same characteristic acceleration (0 in 
cosmic structures of vastly different sizes. (But see 8J.) 

Milgrom has proposed that, instead of dark matter, 
galactic rotation curves signal a modification of gravity at 
very low accelerations j^]. If Newtonian theory predicts 
a gravitational acceleration on, then the gravitational 
acceleration in Milgrom's theory (MOND) is. 
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where 



Vx > 1 
Vx < 1 



(3) 

The MOND interpolating function f{x) is assumed to 
vary smoothly between the stated limits. Many forms 
are consistent with the data. A typical example is. 
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In MOND ao has the status of a new universal constant. 
Its numerical value has been determined by fitting to nine 
well-measured galaxies 



ao 



(1.20 ±.27) X 10-1°^ 
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MOND was designed to explain the TuUy-Fisher rela- 
tion. In the MOND regime of flAr/ao <C 1 one finds, 

y/aoGM 



^ -^vl,= vWGM . (6) 

Making the natural assumption that the total luminos- 
ity L is proportional to the total mass one finds, v'^ — 
aoGM cx L, which is the TuUy-Fisher relation 

MOND not only reproduces the asymptotic rotation 
velocity , it also describes the inner portions of rota- 
tion curves. MOND fits have been worked up for on the 
order of 100 galaxies, using only the mass-to-luminosity 
ratios in stars and in gas as free parameters. The re- 
cent review paper by Sanders and McGaugh sum- 
marizes the data and lists the primary sources. Except 
for a handful of galaxies that show evidence of recent 
disturbance, the fits are excellent. Further, the mass- 
to-luminosity ratios obtained from the fits are in rough 
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agreement with evolution models [ll| . It is especially sig- 
nificant that MOND even fits the rotation curves of low 
surface brightness galaxies 0, 0| , structures for which 
the MOND regime of oat /oq < 1 applies throughout and 
for which no detailed measurements had been made when 
the theory was proposed. Although some dark matter is 
needed to explain the temperature and density profiles 
of large galaxy clusters [lj|, it has been suggested that 
this might be provided by neutrinos without affecting the 
galaxy results . 

The chief problem with MOND is that it does not pro- 
vide a complete theory of gravitation the way that gen- 
eral relativity does. Bekenstein and Milgrom have 
given a satisfactory field theory for the nonrelativistic 
potential whose negative gradient gives the gravitational 
acceleration, a — — V0. If pm is the mass density, the 
Lagrangian for (j) is. 
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where the interpolating function associated with Q 
would be F{x) = Vx + - ln{^/x + ^/l + x). The 
existence of this Lagrangian is very important because 
it establishes that MOND conserves energy, momentum 
and angular momentum. However, this model provides 
no information about the other gravitational potentials, 
or about time evolution. One must therefore make es- 
sentially ad hoc assumptions to test what MOND says 
about lensing 0|, or about cosmology 

What is needed is a generally covariant formulation of 
MOND that includes at least the usual metric. A gener- 
ally coordinate invariant, scalar-metric extension of the 
Bekenstein-Milgrom theory exists but its prediction 
for the deflection of star light is identical to that of gen- 
eral relativity. Without dark matter, this gives far too 
little gravitational lensing from galaxies to be consistent 
with the data |l9j . Scalar- metric models which repro- 
duce MOND and may be consistent with extra-galactic 
lensing data can be constructed, but only at the price of 
allowing the scalar to carry negative energy |20|, |2l| . 

Earlier this year we devised a generally covariant, 
purely metric model of MOND based on a nonlocal effec- 
tive action . Although we were able to reproduce the 
MOND force law, our model also failed to give the extra 
deflection of star light that would be provided in general 
relativity by dark matter. After studying the problem 
we have concluded that it is generic: any stable, gen- 
erally covariant and metric-based theory of gravity that 
reproduces the MOND force law must suffer the same 
problem of too little lensing. The purpose of this letter 
is to give a careful presentation of the argument. We 
do this in section 2. Our model is presented in section 
3 to illustrate the problem. Section 4 identifies the five 
assumptions that lead to the generic problem. We also 
discuss the prospects for abandoning one of them. 



2. The argument: We assume that the gravitational 
force is carried by the metric tensor g^v{x) and that its 
source is the usual stress-energy tensor T^^{x). This im- 
plies that the metric is determined by a set of ten field 
equations having the form. 
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We also assume that the theory of gravitation is generally 
covariant. This implies conservation. 







(9) 



where the semi-colin denotes covariant differentiation. In 
general relativity would be the Einstein tensor, 

G^u, but we make no assumptions about at this 
stage. In particular, we allow it to involve higher deriva- 
tives, and even to be a nonlocal functional of the metric. 

Now recall that deviations between MOND and New- 
tonian gravitation become significant only for very small 
accelerations. These accelerations arise from derivatives 
of the metric. General coordinate invariance allows us to 
choose a coordinate system in which the metric agrees 
with the Minkowski metric ry^j^ at a single point. If 
its gradients are also small — which is the observed 
fact — then the metric can be made numerically quite 
close to rj^i, over a large region. It therefore suffices to 
study equations Q using weak field perturbation theory 
around flat space. That is, we write, 
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and we assume the numerical magnitude of the graviton 
field h^y{x) is small in the MOND regime. 

We wish to consider the weak field expansion of [ry-l- 
h]. In general relativity this begins at linear order, but 
that cannot be the case for any theory which reproduces 
the MOND force law. To see why, consider the gravita- 
tional response to a static distribution of total mass M, 
such as a low surface brightness galaxy, whose density 
is low enough that the equations ^ are everywhere in 
the MOND regime. The MOND explanation © for the 
TuUy-Fisher relation implies that at least one component 
of h^v must scale like \/GM. (For spherical distributions 
it would be the rr component, but this does not mat- 
ter.) How can ||HJ) result in such a dependence? Since the 
right-hand side scales like GM, it follows that at least 
one tensor component of E^,y[rj + h] must go like h^. 

If we assume gravity is absolutely stable then all ten 
components ol E^v[ri + h] cannot begin at quadratic order 
in the weak field expansion. This is because the dynam- 
ical subset of the field equations follow from varying the 
gravitational Hamiltonian. If its variation were quadratic 
then the Hamiltonian would be cubic, and this is not con- 
sistent with stability. We therefore conclude that only a 
subset of the ten components of E^uVl^h] can go like h^ . 
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The desired subset must be distinguished in some gen- 
erally covariant fashion. A symmetric, second rank ten- 
sor field contains two distinguished subsets: its diver- 
gence and its trace. The weak field expansion of the di- 
vergence does vanish a linearized order, but we see from 
conservation ^ that it vanishes to all orders. So the di- 
vergence cannot be responsible for the required term 
and we are left with the trace as the only remaining pos- 
sibility. 



g''''£^,[7^ + h] = 0ih^) 



(11) 



Equation Hll() can explain the TuUy-Fisher relation lO, 
but it means that MOND corrections to general relativ- 
ity can be removed, in the weak field limit, by a local 
conformal rescaling of the metric. 
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To see why, recall that traceless metric field equations 
imply invariance under the transformation (|12]l. (This 
has been known for decades. For a proof, see |23|.) The 
full field equations are not traceless, so neither is the full 
theory conformally invariant. However, because the lin- 
earized field equations are traceless, the linearized theory 
must be conformally invariant. This means that the lin- 
earized field equations only determine the metric up to a 
conformal factor (and, of course, up to a linearized dif- 
feomorphism). The conformal part of the metric is fixed 
by the order term in the trace of the field equations, 
and this is how one component of the weak field h^v con- 
trives to go like y/GM. But these MOND corrections 
to general relativity must be removable by a conformal 
rescaling (|12|l . and possibly a simultaneous coordinate 
transformation. 

This is a disaster for the phenomenology of gravita- 
tional lensing. To see why recall that, for a general metric 
g^^, the Lagrangian of electromagnetism is, 
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where F^^i, = d^Ai, — d^A^. This Lagrangian is invariant 
under a conformal rescaling H12|l . which means that elec- 
tromagnetism is unaware of MOND corrections to general 
relativity in the weak field limit. Hence the deflection of 
star light will be only that predicted by general relativity. 
Without dark matter this gives far too little lensing from 
galaxies |l9j |. That is the generic problem mentioned in 
the title of this paper. 

3. An example: Although the argument we have given 
is completely general, it is illuminating to see how the 
problem arises in a specific model. For this purpose we 
review the model whose failure motivated this study . 
It is based on two nonlocal functionals of the metric. The 
first of these is known as the small potential^ 



if[g] — where 
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We use a spacelike metric with R^y = ^''y^^p — ^'^ p^^^y^ 
r^'p^F'^^^-F^'^^r'^p^. The covariant d'Alembertian is in- 
verted with retarded boundary conditions so that V3[(7](x) 
depends only upon the metric and its derivatives evalu- 
ated on or inside the past light-cone of the point . 

The second nonlocal functional in our construction is 
called, the large potential, 
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To reproduce the MOND force law it turns out that T{x) 
must obey 0|, 



For our model the functional Ep^y [g] is [i^ , 



(16) 
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This form is manifestly covariant. It can be checked that 
it is also conserved. Although £^y[g\{x) is nonlocal in our 
model, it is causal in the sense of depending only upon 
the fields on or inside the past light-cone of the point a;^. 

To prove asymptotic conformal invariance, note first 
that although the weak field expansion of the Ricci ten- 
sor Rpy involves arbitrary powers of the weak fields, it 
begins at linear order. This implies that the weak field 
expansion of the small potential (|14|) also begins at lin- 
ear order. The same applies to the large potential H15|) 
because it is proportional to the small potential in the 
MOND limit oi T{x)^-\x, 



^g] = ^[^PT' 



(18) 



In taking the weak field limit of Spy we can therefore ne- 
glect any products of Rf_iy, ip or <&, such as Gf^y^, 'f'''^,p 
and (p^p(p,y. We can also neglect JF ^ (p^Pip p. Hence the 
weak field limit of Spy is contained in just four terms. 
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Of course these terms involve many higher powers of hpy, 
but they contain all the linear powers. And the terms we 
have kept are exactly traceless, 
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^Dip- R = Q. 

(20) 

Note that tracelessness (and hence conformal invari- 
ance) is not a property of the full theory. In general the 
trace gives, 



g'^'Spy - -6n$ - R[l - 2$] 

+2ip'^^,p + ip^f'ip^pT' ■ 



(21) 
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This goes like h"^ in the weak field expansion. It can 
be shown that the MOND terms in the force law derive 
entirely from the resulting quadratic equation ji^ , in ac- 
cord with the general argument of section 2. 

4- Discussion: The point of making a no-go argument 
is to identify the assumptions which result in the negative 
conclusion. It can then be considered which, if any, of the 
assumptions might be discarded. For the argument given 
above we made the following assumptions: 

• The gravitational force is carried by the metric, and 
its source is the usual stress-energy tensor. 

• The theory of gravitation is generally covariant. 

• The MOND force law is realized in weak field per- 
turbation theory. 

• The theory of gravitation is absolutely stable. 

• Electromagnetism couples conformally to gravity. 

It seems to us that the weakest assumption is absolute 
stability. This is what dictated that only a subset of 
the ten components of £^v[g\ can be quadratic in the 
weak field expansion. If we abandon absolute stability it 
becomes possible that all ten components of £^u[g\ are 
quadratic in the weak fields. 

This may not be as bad as it sounds. It should be un- 
derstood that any theory of MOND necessarily possesses 
two weak field regimes: the ultra- weak limit in which 
MOND applies, and a less- weak regime in which grav- 
ity is weak but MOND corrections to general relativity 
are negligible. It is the latter regime which describes the 
solar system and the interior of our galaxy, so gravity 
would still be stable in these settings. 

The instability could only manifest in regions of very 
low gravitational acceleration. Even there it might be 
self-limiting because the creation of any significant den- 
sity of decay products — whatever they are — would 
likely drive the theory back into the less-weak regime 
in which it is stable. So one might imagine a universe 
that very gradually decays, in the empty regions between 
galaxies, into long wave length particles whose density is 
diffused as the universe expands. If is really a constant 
this decay would only have started recently in cosmolog- 
ical history. And its further progress must be heavily 
suppressed by the intrinsic weakness of the gravitational 
interaction. We therefore conclude it is worth search- 
ing for a generally covariant, metric-based formulation of 
MOND in which all ten components of the field equations 
are quadratic in the weak field limit. 

Even if no viable, generally covariant formulation of 
MOND can be constructed, this would in no way invali- 
date the impressive observational data that has been ac- 
cumulated over many years 0. The absence of an ac- 
ceptable generalization of MOND would mean that this 
data is not explained by an alternate theory of gravita- 
tion, but we wish to stress that the data must still he 



explained. Either the low acceleration regime of gravity 
is ruled by some generalization of MOND or else isolated 
distributions of dark matter evolve towards some hith- 
erto unrecognized attractor solution. Both alternatives 
are fascinating, and we feel the community is greatly in- 
debted to those whose patient labors have drawn atten- 
tion to the problem. 
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